library(raster)
library(sf)
source('https://raw.githubusercontent.com/oharac/src/master/R/common.R')

source(here('common_fxns.R'))

reload <- FALSE

1 Summary

Using the impact category level global maps, aggregate species impacts (for 2011-2013) and intensification (across 2003-2013) to country level. We also aggregate to marine ecoregion level. Values are by count and by proportion of threatened species.

2 Methods

2.1 Set up for calcs

2.1.1 Get IDs for regions

Pull in EEZ raster and MEOW raster, and assemble into a dataframe. Include ocean area for area-weighted averaging.

2.1.2 Define helper functions

Define a function to bring in all rasters (by year) of an impact, and place into dataframe format. From the raster filename, extract the stressor/impact and the year.

2.2 Calculate impact by count and percent

2.2.1 Calculate impacts by Exclusive Economic Zone

Plot the results based on the count and percent of local threatened species impacted by each impact category.

To get the percent by region, divide the mean impacted species by the mean species present. Note that the means can be divided, but the variance of the ratio is not calculated.

2.2.2 Calculate impacts by Marine Ecoregion

Plot the results based on the percent of local threatened species impacted by each impact category.

2.3 Calculate intensification by count and percent

2.3.1 Calculate intensification by Exclusive Economic Zone

Plot the results based on the count and percent of local threatened species under intensification by each impact category.

To get the percent by region, divide the mean intensifying species by the mean species present. Note that the means can be divided, but the variance of the ratio is not calculated.

2.3.2 Calculate intensification by Marine Ecoregion